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We review our knowledge of the phase diagram of QCD as a function of temperature, chemical 
potential and quark masses. The presence of tricritical lines at imaginary chemical potential 
11 = ijT, with known scaling behaviour in their vicinity, puts constraints on this phase diagram, 
especially in the case of two light flavors. We show first results in our project to determine the 
finite-temperature behaviour in the Nf = 2 chiral limit. 
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1. Introduction 

The phase diagram of QCD as a function of temperature T and quark chemical potential jU is 
governed by the interplay of the chiral symmetry and the center symmetry [p. These symmetries 
are exact for zero and infinite quark masses, respectively. Therefore, varying the quark masses 
away from their physical values towards these limits provides useful insight into the behaviour of 
QCD at the physical mass point. 

Similarly, making the chemical potential complex provides enhanced information on the be- 
haviour of QCD at real chemical potential. Since the sign problem which plagues simulations at 
non-zero jx [§] is absent when jU is pure imaginary, the regime of imaginary jU is actually the only 
direction in the complex /I plane where complete, reliable information on the behaviour of QCD 
can be obtained. It turns out that a rich phase diagram as a function of (T,/jL = ijij) and of the quark 
masses (in u = md,m s ) emerges. The critical and tricritical features of this phase diagram, with their 
associated scaling laws, have consequences for the behaviour of QCD at real jU. 

Here, we summarize what is known about this phase diagram and sketch (Fig. 2) a plausible 
scenario, consistent with current numerical simulations augmented with reasonable assumptions 
of continuity of the critical surfaces. We explore in particular the implications for the behaviour 
of QCD in the two-flavor chiral limit (m u = m c i = 0,m s = °°). In that limit, it is widely believed 
that QCD undergoes a finite-temperature, second-order 0(4) chiral transition at /i = 0, which turns 
first-order at a tricritical point for some real /I However, other possibilities exist. At jU = in 
particular, the finite-temperature transition might be first-order. The present numerical evidence is 
inconclusive: using Wilson fermions, 0(4) scaling is preferred [Q], while with staggered fermions 
0(4) scaling has been elusive, and first-order behaviour has also been claimed Note that 
behaviour consistent with 0(4) has been seen with improved staggered fermions, in an Nf = 2 + 1 
setup where the strange quark mass is fixed at its physical value Approaching the chiral limit 
from the imaginary jU direction offers a novel, independent method to help settle the issue. 

2. Three-dimensional Columbia plot 

The thermal behaviour of QCD at jj. = 0, as a function of the quark masses m u =rri4 = m u> d and 
m s is summarized in the well-known Columbia plot Fig. 1 (left). The Nf = 3 chiral symmetry and 
the Z3 center symmetry are achieved in the lower left and upper right corners, respectively. This 
gives rise to first-order transitions. For intermediate quark masses, numerical simulations indicate 
a smooth crossover as a function of temperature. Hence, the first-order regions must be bounded 
by second-order critical lines: the chiral critical line in the lower left corner, and the deconfinement 
critical line in the upper right corner. In the absence of further symmetry, the universality class is 
expected (and has been numerically verified) to be that of the 3d Ising model. The chiral critical 
line joins with the m u ^ = axis at a tricritical point, for a strange quark mass mf 10 which is larger 
than the physical strange quark mass on coarse lattices ^ or smaller when using improved actions 
[0]. The Nf = 2 chiral limit is obtained in the upper left corner. 

When the chemical potential is turned on, the two critical lines sweep critical surfaces as a 
function of }i. For both lines, it has been observed that the first-order region shrinks, as represented 
Fig. 1 (right) [M |8^ ^j. Here, we want to show real and imaginary /I in a single figure. Therefore, 
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Figure 1: (Left) "Columbia plot": schematic phase transition behaviour of Nf = 2 + 1 QCD for different 
choices of quark masses (m u d,m s ) at pi = 0. Two critical lines separate the regions of first-order transitions 
(light or heavy quarks) from the crossover region in the middle, which includes the physical point. {Right) 
Critical surfaces swept by the critical lines as ji is turned on. For light quarks [0, |[ [To|] as well as for 
heavy quarks [^], numerical simulations indicate that the first-order region shrinks as the chemical potential 
is turned on. 



we adopt (jU/T) 2 for the z coordinate: real and imaginary ft appear above and below the /i = 
plane, respectively. Of particular interest is the Roberge-Weiss plane (}i./T) 2 = — (n/3) 2 . We now 
argue that the 3-dimensional phase diagram of Nf = 2 + 1 QCD is likely to be described by Fig. 2. 



3. Phase diagram in the Roberge-Weiss plane 

The two symmetries of the partition function 

Z M =Z(- M ), z(|)=z(|-H^) (3.!, 
imply reflection symmetry in the imaginary /» direction about the "Roberge-Weiss" values /» = 



inT /3(2n + 1) which separate different sectors of the center symmetry [11]. Transitions between 
neighbouring sectors are of first order for high T and analytic crossovers for low T [ 1 1 , [12|, 1 3 ] , as 
indicated Fig. 3 (left). The corresponding first-order transition lines may end with a second-order 
critical point, or with a triple point, branching off into two first-order lines. Which of these two 
possibilities occurs depends on the number of flavors and the quark masses. 

Recent numerical studies have shown that a triple point is found for heavy and light quark 
masses, while for intermediate masses one finds a second-order endpoint. As a function of the 
quark mass, the phase diagram at \i/T = in/3 is as sketched Fig. 3 (middle). This happens for both 
N f = 2 Ql and ^ =3 0. 



If one assumes that the Nf = 2 and Nf = 3 tricritical points are connected to each other in the 
(m Uj d,m s ) quark mass plane, the resulting phase diagram is depicted Fig. 3 (right), with two tricriti- 
cal lines separating regions of first-order and of second-order transitions. This phase diagram is the 
equivalent of the Columbia plot, now at imaginary chemical potential fi/T = in/3. Note that the 
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Figure 2: 3d phase diagram. The vertical axis is (jx/T) 2 , so that real and imaginary chemical potentials 
are above and below the /i = plane, respectively. The "bottom plane" corresponds to the Roberge-Weiss 
transition value ju/T ' = in/3. The thicker red lines are tricritical. Tricritical points marked "2" and "3" have 
been identified for the Nf = 2 [|l4|] and Nf = 3 [|l5[ theories, respectively. The object of the present study is 
the blue line in the "backplane" m s = °° (Nf = 2) joining two tricritical points. 

assumption of continuity of the tricritical lines can be checked directly by numerical simulations, 
since there is no sign problem for imaginary Li. 

Now, as (jLt/r) 2 is varied between zero and the Roberge-Weiss value —(n/3) 2 , the Columbia 
plot must change from Fig. 1 (left) to Fig. 3 (right). Assuming continuity of the critical surfaces 
at imaginary jU, which again can be checked by numerical simulations, the resulting 3-dimensional 
phase diagram is that of Fig. 2. The two red surfaces ("crural" and "deconfinement") are critical. 
They are bounded by lines, among which the following are tricritical: (i) the two lines in the 
(li/T) 2 = — (tt/3) 2 Roberge-Weiss plane; (ii) the line in the m U( i = chiral plane. Note that the 
Nf = 2 (i.e. m s = oo) "backplane" contains two tricritical points on the chiral critical surface: one 
in the Roberge-Weiss plane, the other on the m u c i = vertical axis (see Fig. 2). The location of the 
latter is related to the value of the tricritical strange quark mass. 

4. Tricritical scaling 

In the vicinity of a tricritical point, scaling laws apply. The phase diagram is similar to that of 
a metamagnet, with two external fields: H, which respects the symmetry, and W which breaks it 
(like a staggered and an ordinary magnetic field), depicted Fig. 4 (left). The three surfaces So, S + , S- 
indicate first-order transitions. They meet at a line of triple points L T , depicted by a solid line. They 
are bounded by second-order transition lines, depicted by dotted lines. All four lines meet at the 
tricritical point (T,,H t ). 
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Figure 3: (Left) Generic phase diagram as a function of imaginary chemical potential and temperature. Solid 
lines are first-order Roberge-Weiss transitions. The behaviour along dotted lines depends on the number of 
flavors and the quark masses. (Middle) For Nf = 2 and Nf = 3, the endpoint of the Roberge-Weiss line is 
a triple point (where 3 first-order lines meet) for light or heavy quark masses, and an Ising critical point for 
intermediate quark masses. Thus, two tricritical masses exist. (Right) The simplest assumption is that the 
Nf = 2 and Nf = 3 tricritical points are joined by tricritical lines [Q. This assumption can be checked with 
Nf = 2 + 1 imaginary-pi simulations. 



In our case, the scaling exponents governing the behaviour near the tricritical point are mean- 
field, because QCD becomes 3-dimensional as the correlation length diverges while the temperature 
is fixed to that of the tricritical point, and d = 3 is the upper critical dimension for tricriticality. Of 
course, this implies the presence of potentially large logarithmic corrections to scaling. 

Here, we are interested in the second-order lines S±, corresponding to a departure from the 
symmetry plane = 0. Along these lines, the scaling law is oc \t\ 5 / 2 , where the reduced 
temperature t is measured along the tangent to Lx- For Nf = 2 QCD, tricritical scaling should be 
satisfied near the tricritical points: 

(/) (At/r) 2 = -(7r/3) 2 : thentft~ [{li/T) 2 + (tt/3) 2 ] , t ~ (m M -m^c), so that 

[(H /Tf + (tt/3) 2 ] oc ( mM - ,„ tric )V2 (4 _ 1} 

(«) m u ,d = 0: then // f ~ m u>d , t ~ [{}x/T) 2 - (jU/r) 2 | tric ] , so that 

IBM" [(At/r) 2 -(At/r) 2 | tric ] 5/2 (4-2) 



It is not clear how broad the scaling window is around each of these two tricritical points. The 
two scaling windows might overlap, leading to a very constrained system of equations, with 3 
unknowns (the two constants of proportionality in eqs.(l.l JO) and (jU/^) 2 |tric — "imc having 
been determined already in [|l4|]) and one constraint (continuity of the derivative at the intersection 
of the two scaling curves), leading to a phase diagram as in Fig. 5 (left), which could be determined 
from only two points measured by Monte Carlo. Reason for such optimism can be found in Fig. 4 
(right), where the scaling window around the third Nf = 2 tricritical point, corresponding to heavy 
(u,d) quarks, is shown to extend far into the region of real chemical potential [15, T 
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Figure 4: (Left) Schematic phase diagram of a metamagnet. The external field H T breaks the symme- 
try, while H does not. (Right) For heavy quarks, tricritical scaling in the vicinity of the Roberge-Weiss 
imaginary-jti value extends far into the region of real pL [|l5|]. 



5. Preliminary Nf = 2 results 

Following the above discussion, we performed simulations of Nf = 2 QCD, with staggered 
quarks of masses am q = 0.01 and 0.005 onN t = A lattices, scanning in (/J./T) 2 to determine the 
value of imaginary pt corresponding to a second-order transition. Our observable is the Binder 
cumulant of the quark condensate. Consistent results are obtained from the finite-size scaling of 
the plaquette distribution. The two critical points are shown Fig. 5 (right). Disappointingly, it 
seems impossible to smoothly match two tricritical scaling curves passing through these points. 
Additional masses are needed to determine the critical curve. Nevertheless, assuming convexity 
of the critical curve already constrains the m u ^ = tricritical point to lie at (/J./T) 2 > —0.3. The 
figure illustrates the case where this point lies at pt = 0. It might also lie at (pi/T) 2 > 0, so that the 
/I = chiral transition would be first-order. Additional small-mass measurements are underway 
and will settle this issue. Note that we are simulating two-flavour QCD by taking the square root of 
the staggered determinant, and approaching the chiral limit at fixed, rather coarse lattice spacing. 
This is the wrong order of limits, and is the most likely approach to expose a failure of rooting. 

Finally, our phase diagram Fig. 2 makes it clear that, if the transition in the massless Nf = 2 
theory is 0(4), turning on a real chemical potential (i.e. going up the vertical axis in the back) will 
not make it become first-order. Obtaining such behaviour requires that the chiral critical surface 
bend away from the Nf = 3 chiral point, or that another, non-chiral critical surface appear at large fi. 
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Figure 5: (Left) Strategy used: to determine the critical line, we fix the quark mass and measure the corre- 
sponding critical imaginary chemical potential. The two tricritical scaling lines may match smoothly if the 
scaling windows overlap. (Right) Preliminary results: the scaling windows do not seem to overlap; convexity 
places the m u j = tricritical point at (jx/T) 2 > —0.3. For illustration, it is placed at /X = in the figure. 
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